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Abstract 

On 10 March 2001 the active region NOAA 9368 produced an unusu- 
ally impulsive solar flare in close proximity to the solar limb. This flare 
has previously been studied in great detail, with observations classifying 
it as a type 1 white-light flare with a very hard spectrum in hard X-rays. 
The flare was also associated with a type II radio burst and coronal mass 
ejection. The flare emission characteristics appeared to closely correspond 
with previous instances of seismic emission from acoustically active flares. 
Using standard local helioseismic methods, we identified the seismic sig- 
natures produced by the flare that, to date, is the least energetic (in soft 
X-rays) of the flares known to have generated a detectable acoustic tran- 
sient. Holographic analysis of the flare shows a compact acoustic source 
strongly correlated with the impulsive hard X-ray, visible continuum, and 
radio emission. Time-distance diagrams of the seismic waves emanating 
from the flare region also show faint signatures, mainly in the eastern 
sector of the active region. The strong spatial coincidence between the 
seismic source and the impulsive visible continuum emission reinforces the 
theory that a substantial component of the seismic emission seen is a re- 
sult of sudden heating of the low photosphere associated with the observed 
visible continuum emission. Furthermore, the low-altitude magnetic loop 
structure inferred from potential-fleld extrapolations in the flaring region 
suggests that there is a significant inverse correlation between the seis- 
micity of a flare and the height of the magnetic loops that conduct the 
particle beams from the corona. 



1 Introduction 



Recent developments in the study of flare acoustic emissions Donea and Lind- 
20051 IDonea et cd] [20061 jMoradi et d] [20071 [Martfnez-Oliveros et al] 



sey 



2007 have bolstered the view that seismic emission from flares offers major new 



insights both into flare physics and hehoseismology, ranging from a greatly im- 
proved understanding of flare dynamics and kinematics, to an understanding 
how seismic emission is generated differently by turbulence in magnetic subpho- 
tospheres from inside the quite Sun. 

Sunquakes emanate from compact sources that take only a small fraction of 
the energy released by flares. The surface manifestation of these sources appears 
as circular (or near-circular) waves propagating outward from the solar surface, 
approximately 20-60 minutes after the impulsive phase of the flare. |Donea| 



and Lindsey 120051 considered the possibility that relatively weak flares might 



be able to produce sunquakes and that acoustically active flares may indeed 
be more common than previously thought. This was conflrmed soon after by 
Besliu-Ionescu et al. 2006 following comprehensive helioseismic observations 



of flares using helioseismic holography and the data from Michelson Doppler 
Imager (MDI) onboard the Solar and Heliospheric Observatory (SOHO). 

While the majority of the flare acoustic transients discovered to date have 
been released by the more energetic X-class flares, recently however, a number 
of strong acoustic emissions from M-class flares have been discovered. Donea 



et al. 2006 analyzed the helioseismic properties of the strong seismic tran- 



sient produced by the M9.5 class flare of 9 September 2001, and very recently, 



Marti'nez-Oliveros et al. 2007 performed a comprehensive electromagnetic and 
acoustic analysis of the M7.4-class flare of 14 August 2004 ~ the smallest flare 
known to have produced a detectable acoustic transient prior to the discovery 
of the sunquake reported in this paper. Their findings, along with those of 
Donea and Lindsey 2005 and Moradi et al.\ [2007 , have identified a number of 



distinct observational characteristics that distinguish acoustically active fiares 
from others: 

1. The sites of seismic emission generally coincide spatially with impulsive 
hard X-ray (HXR) and microwave (MW) emissions, suggesting a relation 
to thick-target heating of the chromosphere by energetic particles. 

2. The sites of seismic emission similarly coincide spatially with impulsive 
continuum emission, suggesting acoustic emission associated with extra 
heating and ionization of the low photosphere. 

3. The seismicity of the active region appears to be closely related to the 
heights of the coronal magnetic loops that conduct high-energy particles. 

In this paper, we will examine the last of the above characteristics - intro- 

in greater detail, with further evidence 



duced by Martmez-Oliveros et al. 2007 



from the fiare of 10 March 2001. The evidence we will provide reinforces the 
theory that shorter coronal loops are more likely to be conducive to a more 
rapid injection of trapped, high-energy electrons into the chromosphere at their 
footpoints. This enhances the magnitude and suddenness of the chromospheric 
heating that gives a rise to the intense visible continuum emission seen in all 
acoustically active fiares. This mechanism appears to be a prospective source of 
the energy required to drive a powerful acoustic transient into the solar interior. 
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2 The Helioseismic Signatures 



The helioseismic analysis relies on raw data from SOHO/MDI. The data consist 
of full-disc Doppler, magnetogram, and continuum images in the photospheric 
line Ni i 6768 A, obtained at a cadence of one minute. The Dopplergrams were 
corrected for small effects due to reduced oscillatory amplitudes in magnetic 
regions, following the method outlined in Rajaguru et al. 2006 . We utilize two 



different, but in our opinion, complementary helioseismic techniques to analyse 
the seismicity of the acoustic emission produced by the flare. 

The first method employed was the time-distance technique described by 
Kosovichev and Zharkova^ 11998,. We generate the time-distance plot over a se- 



lected range of azimuths from the primary HXR and magnetic-transient sources, 
in this case -1-135° to -1-225°, in order to gauge the expanding signal from this 
region and compare this signal with a curve that represents the theoretical group 
travel time. The resulting signature, manifested as a "ridge" in the time-distance 
diagram, was significant, but as was expected, appeared to be quite weak (see 
Figure [l}. This is more than likely a consequence of the relatively small energy 
released by the fiare (class M6.7 in X-rays) that produced the sunquake. 

The theoretical curve appears to match the observed ridge with a delay of 
approximately 5 minutes from the time of the fiare maximum. A temporal 



1 Zharkova and Zharkov 


2007 , in our 


According to 


Zharkova and Zharkov| 



case being of slightly longer duration 
|2007j, this delay is due to the time required for the electrons to move along 
the magnetic field lines and hit the upper photosphere or chromosphere. The 
velocity and acceleration of the expanding wave packet was also computed. The 
velocity of the wave-front was calculated to be between 5 and 9 Mm was « 
13 kms~^ and between 29 and 33 Mm, 66.67 kms^^. The mean acceleration 
of the wave-front was also estimated to be « 3.35 kms^^. 

The second method employed in our analysis was computational seismic 
holography, to image the acoustic source of the sunquake. This method has 
been used extensively in the analysis of acoustically active fiares, with great 
success in identifying numerous seismic sources from solar fiares [Donea, Braun 



and Lindsey] p99| [DMlea and Lindsey| |2005l ponea et a?.||2006l 



Moradi et al.\ 



2007| iMartfnez-Oliveros et al. 2007 . Helioseismic holography can be described 
as essentially the phase-coherent reconstruction of acoustic waves observed at 
the solar surface into the solar interior to render stigmatic images of subsurface 
sources that have given rise to the surface disturbance. In general, the acoustic 
reconstruction can be done either forward or backward in time. When it is 
backward in time, we call the extrapolated field the "acoustic egression". In 
the case of subjacent-vantage holography this represents waves emanating from 
surface focus downward into the solar interior that have subsequently refracted 
back to the surface in an annular pupil surrounding the source. For the sake of 
brevity, we direct the reader to Lindsey and Braun 2000 for a more in-depth 



discussion of holographic techniques. 

To assess the seismic emission from the fiare, we computed both the acoustic 
and egression power over the neighbourhood of the active region at one-minute 
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Figure 1: Time -distance plot of the amplitude of the surface ridge averaged 
over curves of constant radius in the azimuth range +135° to +225° are ren- 
dered in gray in both frames. The white curve superimposed in the right frame 
represents the wave time travel for a standard model of the solar interior. The 
time represented as along the vertical axis of the plot is 04:07 UT. 

intervals, mapping them for each minute of observation. It is important to 
distinguish between the "egression" power - wherein each pixel is a coherent 
representation of acoustic waves that have emanated downward from the focus, 
deep beneath the solar surface, and re-emerged into a surrounding annular pupil, 
and the local "acoustic" power - wherein each pixel represents local surface 
motion as viewed from directly above the photosphere. 

The resulting acoustic and egression power movies and "snapshots" (acous- 
tic/ egression power sampled over the solar surface at any definite time) are 
computed over 2 mHz bands, centred at 3 and 6 mHz. The higher frequency 
band has a number of advantages because it avoids the much greater ambient 
noise of the quiet Sun that predominates the 2-4 mHz frequency band and due 
to a shorter wavelength, it also provides us with the images that have a finer 
diffraction limit. 

Acoustic and egression power snapshots at the maximum of the flare are 
shown in Figure [2] In these computations, the pupil was an annulus of radial 
range 15-45 Mm centred on the focus. To improve the statistics, the original 
egression power snapshots are smeared by convolution with a Gaussian with 
a 1/e-half- width of 3 Mm. The egression power images and the continuum 
image are also normalized to unity at respective mean quiet-Sun values. The 
acoustic signature of the flare - consisting of a bright compact source - is clearly 
visible at 6 mHz in the both acoustic and egression power snapshots at 04:05 UT 
(indicated by the arrows in Figure |2|. At 3 mHz the egression and local acoustic 
power snapshots show a less conspicuous signature than at 6 mHz due to a much 
greater background acoustic power at 3 mHz. 

The temporal profiles of the seismic source, seen in the acoustic/egression 
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Figure 2: Egression and acoustic power snapshots of AR 9348 on 10 March 2001 
integrated over 2.0-4.0 mHz and 5.0-7.0 mHz frequency bands and taken at 
the maximum of the correspondence frequency. Top frames show MDI magne- 
togram of the active region (right) at 04:05 UT and a visible continuum image 
at 04:08 UT(left). Second row shows egression power at 3 mHz (left) and 6 mHz 
at the respective maxima. The bottom row show acoustic power. Times are 
indicated above respective panels, with arrows inserted to indicate the location 
of the seismic source. 



time-series in Figure |3] correspond closely with other compact manifestations 
of the flare including significant white-light (WL) emission with a sudden, im- 



pulsive onset as discussed by ,Li, Ding and Liu, ,2005^ and Uddin et at. 



2004 
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Figure 3: The 3- and 6-mHz egression and acoustic power time series, integrated 
over the neighbourhood of the egression power signatures. The vertical Hues 
represent the beginning (04:00 UT), maximum (04:05 UT), and end (04:07 UT) 
times of the GOES X-ray flare. 



The spatial and temporal features of the seismic source observed also coincides 
closely with the HXR signature reported by ILi, Ding and Liul 12005] , indicat- 



ing that high-energy particles accelerated above the chromosphere contribute to 
the generation of the seismic source. We will discuss their observations in more 
detail in the next section. 



6 



3 Multi- Wavelength Analysis 



The multi- wave length properties of the extremely impulsive WLF of 10 March 



2001, have previously been studied in detail by a number of authors [Liu, Ding 



and Fang 2001 Ding et al. 2003 Uddin et al. 2004 Li, Ding and Liu 2005 ; aU 



emphasizing the impulsiveness of the flare and the strong spatial and temporal 
coincidence of the hard HXR emission with the enhanced continuum emission. 



The observations of Uddin et al. 2004 showed that the flare embodied a 



very hard spectrum in HXR, a type II radio burst, and a coronal mass ejection. 
GOES SXR observations classified it as a M6.7 class, beginning at 04:00 UT, 
reaching its maximum at 04:05 UT, and ending at 04:07 UT. A very important 
characteristic of the flare of 10 March 2001 is its duration, which was approxi- 
mately seven minutes, indicating that the physical processes associated with the 



flare also had a very short duration. Uddin et al. 2004 made a detailed study 



of this flare at different wavelengths and determined that all three main phases 
of the flare could be observed clearly in different temporal profiles in HXR at 
different energy bands (Figure |4]). The precursor phase was observed to occur at 
04:03 UT with a duration of 15 seconds, the impulsive phase between 05:03:15 
and 04:03:40 UT, and the gradual phase after 04:03:40 UT. Also, they calcu- 
lated the column emission measure, the spectral index of the flare signal and the 
temporal variation of the temperature. They found that the emission has a non- 
thermal component before 04:04 UT and thermal component after 04:05 UT. 
From the observed profiles, they concluded that a very fast acceleration of the 
electrons occurs during the impulsive phase. 



Uddin et al. 2004 also emphasized the spatial and temporal correlation of 
the HXR source and the continuum emission. They also commented on the 
uncommon change of magnetic flux they detected, concluding that it indicates 
that the WLF was triggered by a new emerging flux that induces a flux cancel- 
lation. As a result, they conclude that magnetic reconnection occurred in the 
upper atmosphere of the sunspot region; thereby high-energy electrons precipi- 
tate along magnetic field lines and deposit energy at the sunspot region, which 
produce the HXR and continuum enhancement. 

The importance of this particular type of spatial and temporal correlation 
between the different types of multi- wavelength signatures described above, in 
the presence of a seismic source, was first identified and discussed in depth by 



Martfnez-Oliveros et al. 12007]. They identified a significant temporal correla- 



tion between the fluxes at different frequencies and energy bands (for the M7.4 
class flare of 14 August 2004) which were seen to be directly related to two elec- 
tron populations - one trapped in the magnetic field, and another precipitating 
into the chromosphere. The highly impulsive character of this flare indicates 
that the trapped population of electrons in the magnetic field was injected into 
the chromosphere very fast. The electrons had no time to thermalize in the 
coronal loop, but were evacuated by rapid precipitation, therefore they did not 
produce a significant emission in MW. Indeed, this type of emission is absent 



in the MW profile reported by Uddin et al. 2004 . The radio emission does 



not show a long exponential decay, implying that high-energy electrons that are 
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Figure 4: HXR and MW time profiles. The HXR fluxes were taken by Yohkoh 
at the L(14-23 keV), Ml(23-33 keV), M2(33-53 keV), and H(53-93 keV) 
channels. The NoRP flux plotted correspond to the 17 GHz channel. 



generally trapped for a signiflcant amount of time in long coronal loops that 
extend to great heights, are evacuated by rapid precipitation in short, low-lying 
loops. 

Li, Ding and Liu 2005 also observed the WL properties of the 10 March 
2001 flare, detecting an infrared continuum enhancement of 4-6% compared to 
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Figure 5: HXR contours of the flare at 04:03:38 UT overlaid over a) MDI 
intensity continuum, b) MDI niagnetogram, c) MDI magnetograni difference 
at the flare maximum, d) Doppler difference at the maximum, e) Ha and f) 
SOHO/EIT at 171 A. The background images all correspond to the same time 
(04:04:01.61 UT). The HXR contour levels are 20, 40, 60, 80, and 90% of the 
maximum emission in the M2 (22-53 keV) channel. The MDI niagnetogram 
neutral line (red line) is overlaid in the frames a), b), c) and d). The blue and 
yellow circles in all frames represent the relative position of the main magnetic 
transients. The seismic source coincides spatially with the blue circle, where 
there is also the HXR emission. 
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pre-flare values. The study of the continuum images shows that the WL source 
is located over the magnetic neutral line and that the source is most likely com- 
posed of the two footpoints of the magnetic loop, which are too close together 
to be resolved by the RHESSI HXR observations. They also detected a HXR 
source near the sunspot. The authors also concluded from their observations 
that the temporal and spatial coincidence of the HXR emission with the con- 
tinuum emission indicates that electron precipitation may have been the main 
energy source of the chromospheric heating, producing the excess continuum 
emission. Furthermore, they suggest that the electron-beam bombardment, 
coupled with radiative back-warming effects, plays the main role in the heating 
of the sunspot atmosphere. This is significant because all instances of seismic 
emissions to date have exhibited very similar WLF characteristics - character- 
ized in particular by the sudden appearance of the WL signature during the 
impulsive phase of the acoustically active flare. 
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Figure 6: Vector magnetogram of the active region taken by the Mitaka Obser- 
vatory (NAOJ) at 00:10:16 UT. The sunquake region is pointed with an arrow. 

The images in Figure |5] show a number of the multi-wavelength signatures 
emitted by the 10 March 2001 solar flare. Frames [Sji and[5]D show the position 
of the magnetic transients, represented by the yellow and green circles, over the 
MDI-intensity continuum and magnetogram respectively. The magnetic neutral 
line is over-plotted (red line) in all frames for reference. FigurejsJ: shows the mag- 
netic difference maps at the time of the maximum of the flare (04:04:01.61 UT). 
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We can clearly see one transient coincides well with the region of HXR emission 
(denoted by the contours), lying across the magnetic neutral lines. In Figure [sji 
we have plotted the Doppler differences for the same time. Here we can see two 
photospheric signatures (spatially coinciding with the magnetic transients) that 
can be associated with surface perturbations of the solar photosphere. We also 
note that observations by Li, Ding and Liu 2005 show that the WL signature 
is composed of two sources, both of them being well correlated (spatially) with 
the magnetic transients. One strong and extended source lies in the region of 
the HXR and seismic source; the second one appears to correlate well with the 
second magnetic transient. 



Uddin et at. 2004 extensively analyzed the temporal and spatial behavior 



of the solar flare. The maximum time in both HXR and MW emission reported 
by them as well as by Li, Ding and Liu 2005] (who undertook a very similar 
analysis) , coincides very well with the maximum of the seismic emission (follow- 
ing the already well known delay of approximately three - four minutes [Moradi 
et al.\ p007| ). Uddin et al. 2004 also discuss the spatial correlation between 



the different sources, showing that all three forms of emissions (WL, HXR and 



MW) are located in the region of maximum magnetic shearing. Chandra et al. 
[2006] , in a similar work, reported the locations of two Ha kernels in the flaring 
region. One of these kernels is (spatially) well correlated with the HXR source 
(observed by Yohkoh) and the observed seismic source, suggesting the precipi- 
tation of electrons in the chromosphere. The second Ha kernel is however not 
correlated with any HXR source, possibly indicating proton precipitation in this 
region (see Zharkova and Gordovskyy 2004 for a discussion about the partial 
separation of electrons and protons into the loop legs). 



4 The Magnetic Field 



The magnetic fleld topology of the active region has also been studied by other 
authors Uddin et al. 2004[ |Li, Ding and Liu 2005 and was correlated with 
other emissions produced by the flare. Using vector magnetograms from the 
Mitaka Solar Observatory (Figure [6|, we can see that the shearing of the mag- 
netic field lines is close to 80° at the location of the seismic source (see the white 
arrow) - which would imply that a vast amount of energy was stored in the mag- 
netic field prior to the fiare. The area where the shearing is significant is very 
small. The seismic source itself is proved to be of a small size of 19 x 25 Mm. 
The magnetic energy released by the flare is used to accelerate particles, heat 
the chromosphere, and also drive the coronal mass ejection [see Uddin et al.\ 



2004 Li, Ding and Liu[ 2005 , and produce the compact seismic source 



In order to verify the magnetic-field configuration of the active region (par- 
ticularly in the corona), we computed the non-linear force free field (NLFFF) 
coronal magnetic field extrapolations of the active region using vector mag- 
netograms from the Mitaka Solar Observatory. The resulting extrapolations 
(seen in Figure [?]) clearly show high-altitude magnetic field lines connecting the 
two leading sunspots of the group, while between the leading and the following 
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sunspots, only low-lying loops are visible (see arrow in Figure [7|. A cornparison 
between the extrapolations with SOHO/EIT images at 17lA (Figure pf) and 
Figure |7| shows that our derived coronal magnetic-field extrapolations are in 
agreement with the observed magnetic field. Because of the close proximity 
of the sunspot to the solar limb and other observational constraints, it is not 
entirely possible to fully reconstruct the complete configuration of the magnetic 
field (in the fiaring region). But nonetheless, we can qualitatively infer the 
overall structure of the coronal magnetic field from our estimates. 



In a closely related work, Chandra et al. 2006 conducted a detailed study of 



the dynamics of 10 March 2001 fiare. As mentioned previously, they identified 
two Ha kernels, with only one kernel (Kl) found to be spatially correlated with 
the HXR emission (see Figure |5^) and therefore with the seismic source. The 
second Ha kernel, labelled K2, has an elongated structure. No HXR emission 
has been correlated with this source, we have also not detected any seismic 
source from this region despite the WL signature present at « 04:04 UT. These 
findings, along with observations of the flaring region made by the SXR telescope 
onboard Yohkoh, led the authors to propose a possible configuration of the 
magnetic field composed of two magnetic loops sharing one footpoint ( "three- 
legged" configuration), and associated with the single HXR source observed by 
Yohkoh. One of the loops appears to be connecting the shared footpoint with an 



opposite-polarity region associated by Chandra et al. 2006 with a secondary. 



stronger, yet distant MW source. The second loop is a low-lying loop connecting 
the shared footpoint with another located inside the region with a high degree 
of magnetic shearing. Furthermore, it is important to state that the two kernels 



observed by Chandra et al. 2006 spatially coincide remarkably well with the 
magnetic transients observed in Figure [5^), with only one of them also being 
well correlated with the HXR emission and the Doppler signature. 

The existence of a relationship between the height of the coronal mag- 
netic loops and the seismicity of active regions has previously been proposed 



in Martinez-Oliveros et al. 2007 . The idea behind this assertion was that 
electrons in short, low-altitude magnetic loops precipitate more effectively than 
long, high-altitude loops because of enhanced scattering by thermal electrons 
ablated from the chromosphere. Electrons whose pitch angles are greater than 
the loss-cone threshold, are trapped in the corona until they are scattered into 
the loss cone. Eventually, these electrons precipitate into the chromosphere and, 
depending on their energy, into the photosphere, transferring efficiently energy 
and momentum to the system. This scattering rate is greatly increased when 
the population of thermal electrons in the loop is large. This generally depends 
on the ablation of chromospheric gas into the corona by the fraction of electrons 
that were initially injected into the loss cone. The volume of a short, low-lying 
loop is much smaller than that of a long high-altitude loop. The electron den- 
sity that results from a given mass of the chromosphere having been ablated is 
thus inversely proportional to this volume. Hence, given these understandings, 
we propose that short, low- lying loops become efficient scattering environments 
promptly greatly expediting precipitation on time scales conducive to seismic 
emission. 
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Figure 7: NLFFF magnetic field extrapolation. The arrow shows the low-lying 
magnetic field region associated with the seismic emission. 



The collapse or relaxation of a high-altitude loop into a low-altitude one due 
to reconnection can greatly expand the loss cone, which would then enhance 



the precipitation distribution if pitch angles were left unchanged [Aschwanden 



2004 . As we understand it, such a collapse facilitates electrons, initially trapped 
in the coronal magnetic field, to precipitate into the chromosphere and photo- 
sphere. Observation in Ha Uddin et al. 2004 of this flare, show the evolution 



of the filaments in the flaring region, changing from a potential configuration 
to a sigmoidal structure due of the high shearing of the magnetic field, with 
a post-flare relaxation of the magnetic field lines also observed in Ha. This 
suggests that the above scenario of electron injection could very well take place, 
making the electron precipitation process much more efficient. 



5 Discussions 

The standard flare scenario divides the flare process into a number of phases. In 
this scenario, the flare particles are accelerated to relativistic or super-relativistic 
velocities in the corona and injected into magnetic field loops whose footpoints 
are in active-region chromospheres. Inevitably, some particles are going to be 
trapped in the coronal magnetic field, while others, those in the magnetic loss 
cone, will precipitate directly into the chromosphere. Eventually the majority 
of the trapped particles are either scattered into the loss cone and precipitated, 
or thermalized (or both) by thermal plasma in the magnetic loop. In the case 
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of the very sudden and impulsive flare of 10 March 2001, the hypothesis is that 
acceleration and injection of particles into the magnetic loop occurred in a short 



period of time Uddin et at. 2004 



This kind of phenomenon can be described using the trapping and injection 

In this model the rate of precipitation 



model proposed by Aschwanden 2004 



of charged particles into the chromosphere is controlled by the relaxation time of 
the system. The aperture angle of the loss cone changes with time, significantly 
opening as the magnetic field collapses to a more potential configuration. It is 
important to note that in this model, the time of acceleration and injection of 
the particles into the magnetic field are almost the same and relatively short 
compared with the precipitation and trapping time. It is fair to assume that if 
the relaxation time is short, the aperture of the loss cone also will change rapidly, 
allowing more particles to reach the chromosphere in a short period of time. This 
depends on eflicient scattering of high-energy electrons into the expanded loss 
cone, which is greatly enhanced by chromospheric ablation of thermal plasma 
into short, low-lying loops. Rapid evacuation of trapped electrons is suggested 
by observations of a rapid decay in non-thermal MW emission. As a general 
rule, thermalization of particles in a magnetic trap is small compared to losses 
due to precipitation. Hence, high-energy electrons evacuated from the coronal 
loop in this way contribute to HXR bremsstrahlung emission substantially as 
well as their counterparts that were initially injected into the loss cone. 

A much more complex model of particle precipitation, in which processes 
such as non-thermal excitation and ionization of hydrogen atoms, and non- 



thermal plasma heating (coulomb and ohmic) is explored by Aboudarham and 



Henoux 1986 , Zharkova and Kobylinskii [1993 , and Zharkova and Zhaxkov 



[2007 . Interestingly, the later show the ohmic heating of the corona by the 
electron beams is so effective that the corresponding particle-induced downward 
propagating shocks are almost depleted of energy, leaving very little energy to 
reach the photosphere and induce any kind of seismic activity. Perhaps this is 
the explanation for why we did not see any seismic sources at the location of the 
Ha K2 kernel in Figure [3j However, we also want emphasize here the possibility 
that photospheric heating also contributes to flare acoustic emission. 

As the whole flaring process occurred relatively rapidly (and given the highly- 
impulsive properties of the 10 March 2001 flare, it is not unreasonable to assume 
so), the solar chromosphere was heated quite suddenly. As we have already seen, 
the multi-wavelength emissions of the flare indicate this; furthermore the strong 
spatial and temporal correspondence between the different types of emissions 
point to radiative back-warming playing a significant role in the heating mech- 



anism. This conclusion was in fact drawn by both Li, Ding and Liu 2005 



and Ding et al. 2003 to explain the origin of the continuum feature of the 
10 March 2001 flare in terms of an "electron-beam-heated flare model", with 
chromospheric radiative back-warming suspected being the chief heating agent, 
originating in the temperature-minimum region. 



The above conclusions, when viewed in conjunction with those of Donea and 



Lindsey 2005 , Donea et al. 2006 , Moradi et al. 2007 , and Martmez-Oliveros 



et al. 2007 , provide direct evidence of flare acoustic emission being driven. 
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in part, by heating of the low photosphere. The basic principle here is that 
the chromospheric radiation further heats up the photosphere, with the result 
being of an optically-thick bound-free absorption, which then introduces a 
pressure transient directly to the underlying medium. The photospheric heat- 
ing hypothesis is well supported by our observations and previous ones - which 
all indicate that instances of flare seismic emission have been characterized by 
a close spatial correspondence between the seismic emission and sudden WL 
emission during the impulsive phase of the flare. Radiative fluxes characteris- 
tic of WL emission seen in all acoustically active flares, if emitted downward 
from the chromosphere (as well as upward), are probably sufficient to heat the 
photosphere a few percent within a few seconds of the onset of the incoming 



radiative flux. This process, described by Machado, Emslie and Avrett 1989 
is called "back- warming" . 



According to rough models described by Donea et at. 2006 and Moradi et 



al. 2007 , such heating (if applied suddenly) should cause a pressure transient 
in the heated layer that drives a seismic transient whose energy flux is of the 
order of those estimated for acoustically active flares. The energy invested into 
the seismic transient is in proportion to 

e-{^f (1) 



a fraction of order l^Ic/Ic times the radiative energy suddenly emitted by the 
flare. We therefore expect acoustic emission due to photospheric heating to be 
inefficient in flares whose WL signatures are weak, diffuse, or not very sudden, 
and this is consistent with examples we have encountered to date. 

In conclusion, it must be stated that to date, no one, single mechanism can 
fully explain the mechanics of flare acoustics and their observational signatures 
- because to do so would be a gross over-simplification of the problem. What 
these results have shown is that the study of flare mechanics, as well as he- 
lioseismology, would greatly benefit from the development of detailed models 
of solar flare-induced seismic emission - from the corona, down to the photo- 
sphere, including modeling of realistic active region subphotospheres. Credible 
models would need to include realistic subphotospheric thermal anomalies to 
represent penumbral and perhaps umbral subphotospheres, and realistic pho- 
tospheric magnetic fields extrapolated to depths of a few hundred kilometers 
beneath the photosphere. The later should also include an account for the 
highly inclined magnetic fields that characterize sunspot penumbra - where a 
significant majority of the seismic emissions observed to date have been detected. 

The authors would like to sincerely thank Profs. Paul Cally, Markus As- 
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